5 POKE55,.:POKE56,56:CLR
10 DV=PEEK(186):IFDV<8THENDV=8
15 POKE53280,0:POKE53281,0:PRINT"[147]"
16 POKE53371,0
25 POKE53272,31
30 AD=49152
35 SYSAD:SYSAD+12
40 BX$="[159][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180][157][180]
42 [153]"LOADOPEN(null)ANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDAND(null)":[158]AD[170]9,15
44 [153]BX$"";[163]38)BX$
45 LC[178]0
46 [153]"OPEN(null)ANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDANDNOT"
49 [153]"":[153][163]1)"PRINT     (null)LEFT$VAL CHR$(null)VALVAL(null) ASC(null)ATN(null)LEFT$ LENATN(null)LEN(null)(null)ATN(null)(null)(null)!     "
50 [143] THE GREEN FLASH
55 :
60 [153][163]2)"CONTLEFT$EIGHT OF (null)BSERVER SYSFT ";:L9%[178]8:[141]1040:H[178]Q9
65 [153][163]2)"CONT(null)ENITH DIST. (null)CCULTING EDGE SYS((null)) ";:L9%[178]3:[141]1040:Z0[178]Q9
70 [153][163]2)"CONTATNEROSOL VALXTINCTION SYS(HAZINESS) ";:L9%[178]4:[141]1040:KD[178]Q9
75 [153][163]2)"CONT(null)CCULT SYS(PEEK)CONTOTTOM OR SYS((null))CONTOP OF (null)UN: ":[151]198,.
80 [161]A$:[139]A$[179][177]"T"[175]A$[179][177]"B"[167]80
82 [139]A$[178]"T"[167]TB[178][171]1:[158]AD[170]9,16:[137]85
84 TB[178]1:[158]AD[170]9,16
85 [153]:[141]830
90 DS[178]1:[143] SUN'S DIST IN A.U.
95 RS[178]1:[143] SUN'S RADIUS IN SOLAR RADII
100 TS[178]5700:[143] READUN'S TEMP KELVIN
105 D[178]0:[143] SUN'S DECLINATION IN DEGREES
110 MP[178]1:[143] PLANET'S MASS IN EARTH MASSES
115 PE[178]1:[143] PLANET'S ROTATION PERIOD IN DAYS
120 R0[178]1:[143] PLANET'S RADIUS IN EARTH RADII
125 [153][163]2)"CONT(null)OUR (null)ATITUDE: ";:L9%[178]6:[141]1040:L[178]Q9
130 P[178]1:[143] SEA-LEVEL PRESSURE BARS
135 [153][163]2)"CONTLENURRENT (null)EMPERATURE SYS(ASC(null)): ";:L9%[178]3:[141]1040:CT[178]Q9
137 T0[178]CT[171]32[172]5[173]9
140 OZ[178]3:[143] OZONE LAYER THICKNESS (MM AT STP)
142 [153]:[153][163]5)"RIGHT$S (null)HIS LENORRECT? OPEN(null)/(null)":[151]198,.
143 [161]A$:[139]A$[179][177]"Y"[175]A$[179][177]"N"[167]143
144 [139]A$[178]"N"[167]42
145 HS[178]12:[143] OZONE HT. (MILES ABOVE SEA LEVEL)
148 [153]"LOAD":[153][163]3)"LENATN(null)LEN(null)(null)ATN(null)RIGHT$(null)CHR$... (null)NE (null)OMENT!":[158]AD[170]9,17
150 AA[178][171]1.3:[143] AEROSOL EXPONENT (UNITLESS)
155 AH[178]1:[143] AEROSOL SCALE HEIGHT (MILES)
160 PP(1)[178]78:[143] NITROGEN % BY VOLUME
165 PP(2)[178]22:[143] OXYGEN PERCENT
170 PP(3)[178]0:[143] CARBON DIOXIDE
175 PP(4)[178]0:[143] HYDROGEN
180 PP(5)[178]0:[143] HELIUM
185 PP(6)[178]0:[143] METHANE
190 PP(7)[178]0:[143] AMMONIA
195 PP(8)[178]0:[143] ARGON
200 :
205 [143] CONVERT INPUT DATA TO CGS UNITS
210 DS[178]DS[172]1.496E[170]13:RS[178]RS[172]6.96E[170]10
215 R0[178]R0[172]6.378E[170]8:H[178]H[172]30.48
220 HS[178]HS[172]1.609E[170]5:AH[178]AH[172]1.609E[170]5
225 D[178]D[172]PI[173]180:L[178]L[172]PI[173]180:Z0[178]Z0[172]PI[173]180
230 MP[178]MP[172]5.976E[170]27
235 PE[178]PE[172]24[172]3600
240 T0[178]T0[170]273.1
245 [158]AD[170]9,17
250 [143] CALCUATE AVERAGE MOLECULAR WEIGHT
255 MW[178]0
260 [129]IG[178]1[164]8:MW[178]MW[170]MX(IG)[172]PP(IG)[172].01:[130]IG
265 :
270 [143] CALCUATE ATMOSPHERIC PARAMETERS
275 G[178]6.67E[171]8[172]MP[172](R0[174][171]2)
280 S[178]8.31E[170]7[172]T0[173](G[172]MW)
285 LA[178][171].4[172]T0[173]S
290 :
295 [143] CHOOSE A COLOR
300 [129]IC[178]1[164]11
305 W[178](IC[172]500[170]2500)[172]1E[171]8
310 [143] CALCUATE INDEX OF REFRACTION
315 N0[178]0
320 [129]IG[178]1[164]8
325 NS[178](A(IG)[172].00001)[172](1[170](B(IG)[172]1E[171]11)[173](W[174]2))
330 N0[178]N0[170]NS[172]PP(IG)[172].01[172]P[172]273.1[173]T0
335 [130]IG
340 N[178]1[170]N0[172][189]([171]H[173]S)
345 N0[178]N0[170]1
350 :
355 [143] CAL REFRACTION & AIR MASS FOR EACH COLOR
360 [141]795
365 [139]RF[177]RX[167]RX[178]RF
370 [139]RF[179]RN[167]RN[178]RF
375 R(IC)[178]RF:XR[178]AM
380 [158]AD[170]9,17
385 [143] CALC AIR TRANSMISSION AT OCCULTING OBJECT
390 KR[178]1.13E[171]11[172]((N[171]1)[174]2)[172](W[174][171]4)
395 KA[178]KD[172]((W[173].000055)[174]AA):CZ[178][190](Z0)
400 XA[178]1[173](CZ[170].8[172][186](AH[173]R0)[172][189]([171].38[172]CZ[172][186](R0[173]AH)))
405 K0[178]KS(IC)[172]OZ[173]3
410 X0[178]1[173][186](1[171]([191](Z0)[173](1[170]((HS[171]H)[173]R0)))[174]2)
415 KX[178]KR[172]XR[170]KA[172]XA[170]K0[172]X0:[139]KX[177]90[167]KX[178]90
420 TR(IC)[178]10[174]([171].4[172]KX)
425 :
430 [143] CALCULATE SUN'S FLUX ABOVE ATMOSPHERE
435 BB[178].0000374[172](W[174][171]5)[173]([189](1.43[173](W[172]TS))[171]1)
440 FA(IC)[178]((RS[173]DS)[174]2)[172]BB[172]1E[171]8
445 [130]IC
450 [158]AD[170]9,17
455 [143] PICK RELEVANT TIMR RANGE FOR CALC
460 CZ[178][190](Z0[170]RN[170]TB[172]RS[173]DS)
465 HC[178](CZ[171][191](D)[172][191](L))[173]([190](D)[172][190](L))
470 HA[178][193]([186](1[171]HC[172]HC)[173]HC):[139]HC[179]0[167]HA[178]HA[170]PI
475 T1[178]HA[172]PE[172].5[173]PI
480 CZ[178][190](Z0[170]RX[170]TB[172]RS[173]DS)
485 HC[178](CZ[171][191](D)[172][191](L))[173]([190](D)[172][190](L))
490 HA[178][193]([186](1[171]HC[172]HC)[173]HC):[139]HC[179]0[167]HA[178]HA[170]PI
495 T2[178]HA[172]PE[172].5[173]PI
500 DT[178](T2[171]T1)[173]50
505 T0[178]T1[171]DT
510 [153]"LOADCLR LEN(null)FOR(null)SYS(null)PRINT(null)CONT(null)- WAIT   SYS (null)IME- "
515 [153]"PRINT    -(null)-WAIT        -(null)-WAIT       -(null)-MAG. "
520 [143] MAIN LOOP STARTS HERE
525 [129]IT[178][171]1[164]19
530 T[178]T1[170]IT[172]DT
535 TT[178]T[171]T0
540 :
545 [143] CYCLE THROUGH WAVELENGTHS
550 [129]IC[178]1[164]11
555 [143] CALC ALTITUDE OF THE SOLAR DISK FOR THAT COLOR
560 HA[178]2[172]PI[172]T[173]PE
565 CZ[178][191](D)[172][191](L)[170][190](D)[172][190](L)[172][190](HA)
570 ZA[178][193]([186](1[171]CZ[172]CZ)[173]CZ):[139]CZ[179]0[167]ZA[178]ZA[170]PI
575 ZA[178]ZA[171]R(IC)
580 :
585 [143] CALC THE FRACTION OF DISK VISIBLE
590 CC[178](ZA[171]Z0)[172]TB[172]DS[173]RS
595 [139]CC[177]1[167]CC[178]1
600 TH(IC)[178][193]([186](1[171]CC[172]CC)[173]CC)
605 [139]CC[179]0[167]TH(IC)[178]TH(IC)[170]PI
610 FR[178](TH(IC)[171]CC[172][191](TH(IC)))[173]PI
615 :
620 [143] CALC VISIBLE FLUX FROM SUN
625 F(IC)[178]FA(IC)[172]TR(IC)[172]FR
630 [130]IC
635 :
640 [143] COULD PRINT OUT FLUX F(IC) [ERG/CM2/S/A]
645 [143] AS WELL AS THE SUN'S STELLAR MAGNITUDE
650 [143] M(IC)=-2.5*LOG(F(IC))/B0-(M0(IC)+20)
655 [143] AS A FUNCTION OF THE COLOR IC WHERE
660 [143] 1=3000A, 2=3500A, 3=4000A, ...., 11=8000A
665 :
670 XX[178]0:YY[178]0:ZZ[178]0
675 [129]IC[178]1[164]11
680 XX[178]XX[170]X1(IC)[172]F(IC)
685 YY[178]YY[170]Y1(IC)[172]F(IC)
690 ZZ[178]ZZ[170]Z1(IC)[172]F(IC)
695 [130]IC
700 [139]XX[170]YY[170]ZZ[178]0[167]780
705 X[178]XX[173](XX[170]YY[170]ZZ)
710 Y[178]YY[173](XX[170]YY[170]ZZ)
715 [139]Y[179].01[167]780
720 PH[178][193]((.333[171]X)[173](Y[171].333))
725 [139](Y[171].333)[179]0[167]PH[178]PH[170]PI
730 CO$[178]"CLR (null)RIGHT$(null)(null)VAL(null) "
735 [139]PH[179]4.229[167]CO$[178]"CONT PEEK(null)(null)VAL "
740 [139]PH[179]2.580[167]CO$[178]"PRINT CHR$(null)VALVAL(null) "
745 [139]PH[179]1.106[167]CO$[178]"SYS (null)VAL(null)(null)(null)(null) "
750 [139]PH[179].359[167]CO$[178]"FOR (null)(null)ATN(null)CHR$VAL "
755 [139]PH[179]0.000[167]CO$[178]" (null)VALSTR$ "
760 V[178][171]2.5[172][188](.034[172]YY[173]2.54E[171]10)[173]B0
765 [153]CO$,TT
770 [153]X;Y;V:[158]AD[170]9,18
775 LC[178]LC[170]1:[139]LC[178][177]10[167][141]1100
780 [130]IT
785 [141]3000
790 [153]"LOAD":[137]42
795 [143] REFRACTION SUBROUTINE (QUICK VERSION)
800 N[178]1[170](N0[171]1)[172][189]([171]H[173]S)
805 BS[178]R0[172][191](Z0)[173](R0[170].5[172]S):BE[178][193](BS[173][186](1[171]BS[172]BS))
810 R5[178]N[172]BS:RF[178][193](R5[173][186](1[171]R5[172]R5))[171]BE:CZ[178][190](Z0)
815 AM[178]1[173](CZ[170].8[172][186](S[173]R0)[172][189]([171].38[172]CZ[172][186](R0[173]S)))
820 [142]
825 :
830 [143] CONSTANTS AND DATA STATEMENTS
835 RX[178][171]1:RN[178]1
840 PI[178][255]:B0[178][188](10)
845 [134]A(8),B(8),PP(8),MX(8),KS(11),M0(11),X1(11)
850 [134]Y1(11),Z1(11),R1(11),G1(11),B1(11),M(11)
855 [134]R(11),TR(11),F(11),FA(11),LF(11),TH(11)
860 [129]I[178]1[164]8:[135]A(I):[130]I
865 [131]29.1,26.6,43.9,13.6,3.5,43,37,27.9
870 [129]I[178]1[164]8:[135]B(I):[130]I
875 [131]7.7,5.1,6.4,7.5,2.3,6,12,5.6
880 [129]I[178]1[164]8:[135]MX(I):[130]I
885 [131]28,32,44,2,4,16,17,39.9
890 [129]I[178]1[164]11:[135]KS(I):[130]I
895 [131]3.2,.005,0,.001,.012,.031,.044
900 [131].023,.008,.005,.001
905 [129]I[178]1[164]11:[135]M0(I):[130]I
910 [131]1.19,.95,.52,.46,.74,1.06,1.33
915 [131]1.6,1.88,2.1,2.3
920 [129]I[178]1[164]11:[135]X1(I):[130]I
925 [131]0,0,.014,.336,.005,.433,1.062
930 [131].283,.011,0,0
935 [129]I[178]1[164]11:[135]Y1(I):[130]I
940 [131]0,0,.0004,.038,.323,.995,.631
945 [131].107,.004,.0001,0
950 [129]I[178]1[164]11:[135]Z1(I):[130]I
955 [131]0,0,.068,1.773,.272,.009,0,0,0,0,0
960 [129]I[178]1[164]11:[135]R1(I):[130]I
965 [131]0,0,.025,-.039,-1.16,.097
970 [131].847,.989,1,0,0
975 [129]I[178]1[164]11:[135]G1(I):[130]I
980 [131]0,0,-.011,.0218,1.39,.905,.154
985 [131].001,0,0,0
990 [129]I[178]1[164]11:[135]B1(I):[130]I
995 [131]0,0,.986,1.017,.778,-.0025
1000 [131]-.0012,0,0,0,0
1005 [142]
1040 Q9$[178]"":[151]198,.
1045 [161]A$
1050 [151]646,[187](1)[172]15[170]1:[153]"*CMD";:[139]A$[178]""[167]1045
1055 [139]A$[178][199](13)[167][153]" ":Q9[178][197](Q9$):[142]
1060 [139](A$[178][199](20)[175][195](Q9$))[167]Q9$[178][200](Q9$,[195](Q9$)[171]1):[137]1090
1065 [139][195](Q9$)[177][178]L9%[167][158]AD[170]9,20:[137]1045
1070 [139](A$[177][178]"0"[175]A$[179][178]"9")[176]A$[178]"."[167]1080
1075 [137]1045
1080 Q9$[178]Q9$[170]A$
1085 [153]""A$;:[158]AD[170]9,16:[137]1045
1090 [153]" CMDCMD CMD";:[137]1045
1100 [153][163]6)" ATNNY KEY TO LENONTINUE. ":[151]198,.
1105 [161]A$:[139]A$[178]""[167]1105
1110 LC[178]0
1115 [153]"LOADCLR LEN(null)FOR(null)SYS(null)PRINT(null)CONT(null)- WAIT   SYS (null)IME- "
1120 [153]"PRINT    -(null)-WAIT        -(null)-WAIT       -(null)-MAG. "
1125 [142]
3000 [153]:[153][163]8)"OPEN(1OPEN) (null)RY ANOTHER ONE":[158]AD[170]9,16
3010 [153][163]8)"OPEN(2OPEN) (null)O (null)(null)RIGHT$(null)VAL(null)(null)VAL (null)ENU
3020 POKE198,0
3030 GETA$:IFA$<"1"ORA$>"2"THEN3030
3040 IFA$="1"THENRETURN
3050 SYSAD+15
3060 PRINT"[147][144]LOAD"CHR$(34)"B.UNIVERSE IV"CHR$(34)","DV
3070 PRINT"RUN28"
3080 POKE631,13:POKE632,13:POKE198,2:END
10000 D=PEEK(186):N$="GREEN FLASH":OPEN15,D,15,"S0:"+N$:CLOSE15:SAVEN$,D:END